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Abstract. It is well known that the solutions to the non-viscous Burgers 
equation develop a gradient catastrophe at a critical time provided the initial 
data have a negative derivative in certain points. We consider this equation 
assuming that the particle paths in the medium are governed by a random 
process with a variance which depends in a polynomial way on the velocity. 
Given an initial distribution of the particles which is uniform in space and with 
the initial velocity linearly depending on the position we show both analytically 
and numerically that there exists a threshold effect: if the power in the above 
variance is less than 1, then the noise does not influence the solution behavior, 
in the following sense: the mean of the velocity when we keep the value of 
position fixed goes to infinity outside the origin. If however the power is larger 
or equal 1, then this mean decays to zero as the time tends to a critical value. 



1. Introduction 

The non-viscous Burgers equation is perhaps the simplest equation that models 
the nonlinear phenomena in a force free mass transfer. It has the form 

ut -t- (u, V) u = -/3u, (1.1) 

where u{x,t) — (wi,...,w„) is a vector-function W^^^ M", /3 > is a constant 
friction coefficient. Consider the Cauchy data 

u{x,0) ^uo{x). (1.2) 

Problem (1.1), (1.2) has an implicit solution 

u{t, x) = e~^* uo{x - i (e''* - l)u{t, x)), 
P 

for /5 > and 

u(i, x) = uq{x — tu{t, x)), 

for P = 0. 

In several cases we can obtain an explicit solution. For example, if 

uo{x) — ax, a G K, (1-3) 



one easily gets 



"(^'•^) = TT^(r3-^' /5>0' (1-4) 



Date: August 5, 2008. 

1991 Mathematics Subject Classification. 35R60. 

Key words and phrases. Burgers equation, random position of particle, gradient catastrophe. 
Supported by DFG 436 RUS 113/823/0-1. 

1 



2 



ALBEVERIO, ROZANOVA 



and, for /3 = : 

= (1.5) 

Thus, if a < — /3, the solution develops a singularity at the origin as t — > T, < 
T < oo, where 

T=4ln— for/3>0, T = --, for /3 = 0, a ^ 0. (1.6) 

p a + p a 

This phenomenon is called the gradient catastrophe. It is well known (see [l])that a 
viscous perturbation of form ctAm, a > 0, entails a globally in time smooth solution 
(at least for bounded initial data). An exceptional case is exactly given by a solution 
which is linear in x as mentioned above, which does not feel the viscous term. 

Our main question is: can a stochastic perturbation suppress the appearance of 
unbounded gradients? 

We can introduce the Lagrangian coordinate x{t) to label a point which moves 
together with the medium, that is ^^^^ — u{t,x{t)) :— ui{t). Thus, x = x(t) is 
the equation for the particle path, when the particle moves along the Burgers fluid. 
Equation (1.1) is equivalent to the following system of ODE: 

i;(i) = ui(0, ui(t) = -/3ui(t) (1.7) 

Further on we will omit the index 1. 

In the theory of stochastic dynamical systems one often considers a stochastic 
perturbation of the velocity, which leads to the appearance of a white noise in the 
second of equations (1.7). The problem of solving such equations was investigated 
in many works (see, e.g. [2], [3], [3], [S], p]). This type of stochastic perturbation 
corresponds to the stochastically forced Burgers equation, or in the language of 
physicists, Burgers turbulence. This has been an area of intensive research activity 
in the last decade (see e.g. [7], and for a very recent review [8j, and references 
therein). 

The behavior of the gradient of velocity was studied earlier in other contexts 
in [S], |10| . but this problem is quite different from the problem considered in this 
paper. 

In the present paper we consider a medium with random particles paths, more 
precisely, described by a 2 x n dimensional Ito stochastic differential system of 
equations 

dXk{t) = Uk{t) dt + a\U{t)\P diWk)u 

dUk{t) ^ ~l3Uk{t)dt, fc = l,..,n, (1.8) 
X{0) = X, C/(0) = u, t>0, 
where {X(t), U (t)) runs in the phase space R" x R", a > and p > are constants, 
{W)t — {W)k,t, k ^ 1, ...,71, is the n - dimensional Brownian motion. We remark 
that for p < 1 for initial distributions of x and u from the class L^(R") one can 
guarantee a global existence of a unique solution to (1.8)pT]. 

Let us denote by u{t, x) the mean of the velocity U{t) at time t when we keep 
the value of X{t) at time t fixed but allow U{t) to take any value it wants fe.g.|12|) 
and chose this function for comparison with the solution to the solution of the 
non-viscous Burgers equation. 

We can interpret system (1.8) also as follows: assume that we measure the 
position of a particle with an error depending on its velocity and then try to restore 
the velocity. If the coefficient p increases, the error for large velocities increases, too. 
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It is natural to expect that the mean of difference between two "very indefinite" 
neighbor coordinates necessary to calculate the velocity tends to zero. 

Can we hope to extract from our measurement a realistic information on such 
critical phenomena as the blow up occurring in a medium described by the Burgers 
equation associated with (1.8)? As we will see, at least for uniform initial distri- 
bution of particles (in the sense of Sec. 3, see below) the answer depends on the 
exponent p, namely, if p > 1, the information gets lost. The threshold value p = 1 
is not astonishing. The sub- linear rase of the drift and diffusion (p < 1) coefficients 
warrants the global existence of the SPDE solution provided the initial distribu- 
tions of the particle positions and the velocities are square integrable. Therefore 
one can hope that the solutions to the SPDE with a sub-linear diffusion coefficient 
behave in a some "predictable" way. 

For example, for /3 = the function u{t, x) demonstrates the following behaviour 
near the critical time T. For p e [0, 1) 

u{t,x) = " ^ X + o( "^t j, t^T, xeM"; 

ioY p > 1 

u{t,x) = -C\xf-^ x(^l- + o(l-|), t^T, xy^O, 

where C is the positive constant depending in particular on the dimension of space, 
however, 

u{t,x) = t X + o{\x\) , t€[0,T),x^O. 

1 — y 

In other words, for p G [0, 1) we see that u{t,x) ~ u{t,x), t ^ T at every point 
X G R", where u{t,x) is given by (1.5), i.e. u{t,x) keeps the property of solutions 
to the inviscid Burgers to have a gradient catastrophe. For p > I outside the origin 
a; = the function u{t, x) ^ 0, t ^ T, and the features of the solution of the 
non-perturbed equation fail. A jump is being formed at the origin x = but the 
height of this jump at the same time tends to zero as t ^ T. 

The paper is organized as follows. In Sec. 2 we find an exact solution to the 
Fokker-Planck equation with special initial data and derive an integral formula for 
u. In Sec. 3 we consider the simplest case p = 0, where for several initial distribution 
it is possible to find u exactly. In Sec. 4 we formulate the main theorem for the case 
of uniform initial distribution concerning asymptotics of u near the critical time T 
and the origin x = and prove it. Here we present also the results of computations 
given directly according to the integral formula for u that illustrate the fidelity of 
our asymptotic formulas. In Sec. 5 we present certain results, both asymptotic and 
numeric, concerning a Gaussian initial distribution and show that for < p < 1 the 
behaviour of u is close to the behaviour of the same function at p = 0, where the 
exact formula can be derived. In Sec. 6 we discuss the question on the "observable" 
density and "induced" velocity. In Conclusion we sum up our results and discuss 
the question on a possibility to construct a PDE such that u is its solution. We 
also argue possible applications of the results. 
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2. Exact solution to the Fokker-Planck equation for special initial 

DATA 

The Fokker-Planck equation associated to (1.8) for the probabihty density in 
position and velocity space P — P(t, x. ii) has the form 



dP{t, X, u) 



:i2 



- \u 



2U,|2p, 



fr[\duk J 2 dxl 
subject to the initial data 



P{t,x,u), 
(2.1) 



P{0,x,u) = Po{x,u). 
Let us set := L > 0. Thus, 

/ uP{t,x,u)du 

u{t,x) = lim ^\ t>0,xGflL, (2.2) 

L-yoo J P{t, X, U) du 

provided the limit exists. 
If we choose 

n 

Po{x, u) = 5{u - uo{x)) fix) = n Huk - iuo{x))k) fix), (2.3) 

;c=i 

with an arbitrary sufficiently regular f{x), then 

u{0, x) = uo{x). 

The function fix) has the meaning of a probability density of the particle positions 

in the space at the initial moment of time. 

Moreover, formula (2.2) can be investigated for functions fix) which are not 
necessarily probabilities densities. 

Let us choose 

uoix) = ax, a <Q (2.4) 

as initial data of the non-perturbed Burgers equation. One can see from (1.4), 
(1.5) that the gradient of the solution become unbounded as i ^ T. Thus, we are 
interested in the behavior of w(t, a;) comparing with the solution uit,x) to (1.1). 

We apply formally the Fourier transform in the variable x to (2.1), (2.3) (2.4) to 
obtain for P = Pit, A, u) 

dP dP -2 



= ^3 E + (/5 - Y I - * (A, u)) p, (2.5) 

P(0,A,w) = ^_e-^^/f-V (2.6) 

Equation (2.5) is of the first order, therefore the Cauchy problem (2.5), (2.6) for 
the function P(t, A, u) can easily be solved. Thus, for /3 = 

P(^t,X,u) = ^^""'^ ^-#|.r''|Art-i(A,.)(^+t) (2.7) 

^ ' i\a\^)^ ^ ' 

for /3 > 

Pit,\u) = ^ ("^^'/") c/3t-^l»l^''|Ap(e^'-'-l)-i(A.u)((^ + ^)e/^'-i)_ ^3.8) 

(|a|-\/27r)" 
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Further, the inverse Fourier transform gives for f3 — 

f(u/a) iu(l;+t)-^\^ 

P(t,x,u) = JAW e" '-'^M'" , t>0, (2.9) 

and for /J > 

P(t,x,u) = ; ^ ' ' = e 2<,-i„i-p (.-/^pt-i) ,f>0. 

(|a|c7V27r(e2/5p* - l)//3i)" |u|f" 

(2.10) 

It is easy to see that the limit for /3 — > of (2.10) gives (2.9). 

Now we substitute P{t,x,u) in (2.2) to get an integral representation of the 
quantity expectation u{t,x). It is important to note that for every even function 
f{x) the upper integral in (2.2) vanishes at x = or t = T, hence u{t, 0) = u{T, 0) = 
0. 

However, representation (2.2), (2.9) (or (2.10)) does not allow to see the be- 
haviour of u. Therefore we firstly consider the situations where the integrals in 
(2.2) can be explicitly computed (Sec. 3) or analyze the asymptotics in certain im- 
portant points (Sec. 4). 

3. Case of constant noise variance 

Firstly wc set p = 0. In this section we concentrate on a special choice of the 
initial probability density f{x), allowing to obtain an explicit formula for u{t,x) if 
the initial data have the form (2.4). 

Let us consider uniform initial distribution of particles as follow: 

f{x) = / = W ' ^ f (3.1) 
( 0, otherwise. 

Applying (2.2), (2.9), (2.10) we readily calculate 

U{t,x)^\ l + fll-e-f*)' ^<^. (32) 

[0, t = T 



for /3 > and 



t < T 



u{t,x) = { ^ (3-3) 



for /3 = 0. We also see that (3.3) results in the limit /3 ^ from (3.2). 

It is enough to compare these expressions with (1.4), (1.5) to see that the white 
noise with a constant variance does not influence the mean of the velocity. 
Thus, at any point (T, a;), x ^ 0, the function u{t, x) has a discontinuity (as for the 
case without noise). 

We can also compute the variance 

J {u — uY P{t, X, u) du 
v{t,x) := lim , i>0, xeM". (3.4) 

L^oo J P{t,X,U)du 

We have 

. , ^ at - _ 
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Thus, as f — > T, the denominator of v{t, x) becomes zero and the possible values of 

U (t) at any x dissipate over the space. 

A behavior of the mean of velocity which contrasts with the previous one can be 

-^j exp{—k^x^), A: > 0. It is easy to compute that 

in this case for /3 = one has 

\ (1 + Oit)ax ,„ 

-u \ 
v{t, X) 



The denominator does not vanish for any fixed x, a, a, so u{t,x) 0, v{t,x) — > 
0, t — > oo, and u{t,x) is continuous at any point (T,x). The minimal value of 

divu(t, x) is attained at the time i* = ^ (^^^^cr^ — 1 j < T. Moreover, > only 
if ak < 



Computations made for some special classes of f{x) allow to suggest that a simi- 
lar behavior is provided by the function u{t, x) if f{x) = f{\x\) and / \x\^ f{\x\) dx < 

OO. For example, for the class of f{\x\) = jf^^^sy, n = 1 we have the analytical 
result 

uit,x) ^ —- — — —prxil + at), as a; ^ 0, s e N, s > 2. 

(2(s — 1) — IjK^ 

A numerical study suggests that for n = 1 as a; — > 

fors>l {{(t, x) const a; (1 + ai), 

forse[^,l] a;) consta;, 

1 const 

tor s < - u(t, x) ~ X . 

2 ^ ' ' l + at 

Here the constants do not depend on t. 



4. Asymptotic behavior for uniform initial distribution 

We consider again f{x) given in (3.1). If p > 0, formula (2.2) does not allow to 

compiitc u{t, x) explicitly. Thus, we need to extract from this formula an informa- 
tion that allows us to answer the main question of this article: does the stochastic 
perturbation suppress the singularity formation? 

We will show that there exists a critical value of the parameter, p = 1, such that 
for p < 1 the mean of the velocity behaves very closely to the velocity in the case 
p = (as in the non stochastically perturbed case for t < T). In contrast, for p > 1, 
u{t, x) vanishes a&t^T. 

We deal with the case /? = to avoid cumbersome formulas, the results for /? > 
will be qualitatively the same, we present them in Remark 4.1. 
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Theorem 4.1. The mean u(t,x) of the random variable U(t), given the position 
X{t), where X{t), U{t) solve the SDE (1.8), provided u{Q, x) = ax, a <Q, is given 
by formula (2.2), where P{t,x,u) is given by (2.9), (2.10). 

If (3 = and initially the particles are distributed uniformly in the sense (3.1), 
the asymptotic behaviour of u{t, x) for t ^ T can be analyzed explicitly. 

Namely, for any p € [0,1) the mean u{t,x), being equal to zero at any point 
X G W^, t = T, is discontinuous at every such point if x ^ 0. More precisely, 

• forp = the mean u{t,x) coincides with the solution to the problem (1.1), 
(2.4) for t < T; 

• for p € (0, 1) the asymptotics 

u{t, x) = ^ X + o( ^ t ) ' t^T, xeM."', 

1— J? \1 — j; J 

takes place. 

• For p > 1 at any x G M", x ^ 0, \u{t,x)\ — > as i — > T, however, 
divxu{t, x) — > oo, a; = 0, f — > T. More precisely, for any x ^ 0, 

u{t,x) = -C\x\'-^x(^l-^^ + o{l-^), t^T. 

where C is the positive constant given in (4-4) > and for any t € [0,T) 

cx 

u{t, x) = J- X + o {\x\) , a; — > 0. 

^ ~ T 



Proof. Let us set e := 1 — ^ t, where the critical time T is introduced in (1.6) 
(e e (0, 1]). We will write below t{s) = (1 - e)T. 



I 



Firstly we compute asymptotics of (2.2) near the critical time t = T {e = 0). 
Let us fix a; ^ 0. We have the following expansion as e — > 0,p > 1 

\aL\ 



uP{t{e),x,u)du = f^'f^U-^l^ j eW^^^^-^I__^dHxe+o(e). 

(4.1) 



/ 



/ I \ " 2 

P{t{e),x,u)du = /^a;„ 1^-^ / |„ J-n+i d\u\+0{e), 



(4.2) 

where LOn is the area of surface of (n — 1) - dimensional sphere, n > 2; wi = 2. 
Recall that for fixed x these integrals are functions of e only. 

Both integrals in (4.1) and (4.2) converge as L ^ oo and can be expressed 
through Gamma-functions, so according to (2.2) we get 



u{t,x) = 



£X\X 
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'""^^^^(\/^)" (7^)^""^^"^ r((n + 2)(p-l)/2p) +oie) 

(v^)""^ r(n(p-l)/2p) +0(.) 

2(l-t>) 

= -C£a;|a;|^^ +o(e), e ^ 0, (4.3) 
with the positive constant 

1 / r^"^ r((n + 2)(p-l)/2p) 

r(„(.-i)/..) ■ 

If p < 1, the integrals in (4.1) and (4.2) diverge as i — > oo. Now we take into 
account that 

/ P{t{e),x,u)du= h I ^ Ffe(a;,|M|)£'=d|M|, 



ffc(a;, |«|) = e^^W ^ 



s=0 ' ' 



|aL| 

/ ^ oo 
uP{t{e),x,u) du = ex / ^ Gfc(a;, Im]) e*^ rf| 
Oi ^=0 



k 7. i.|2s 



g,(.jnl) = e-^^ J^^;J^i^'_„_^ . 



s=0 

with some constant bks and Cfeg- Further, 



i(f(e)),a;) = ex hm 



/c=0 



fe=0 



= a;£ Urn — ■ = (4.5) 

/ e^j^rfC(l + 0(£)) 





^+0(£), p=l, 

— OO sgn.T, p < 1. ' 
where we use de L'Hopital's rule to compute the ratio of the divergent integrals in 
(4.5). 

Thus, for p = 1 wc get (4.3) again. 

For p < 1 we are going to obtain a more specified result. Recall that 

P{t(e\x,u) = (4.6) 
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For the third factor as e — > we have: 

|2 



The second factor in (4.6) secures the convergence of integrals. Calculations of both 
integrals in (2.2) with the use of the Maple environment allow to obtain explicit 
formulas for rational p, the result being expressed through special functions (Bessel, 
Gamma and hypergeometric functions). We do not quote here this formula as it is 
very cumbersome. The simplest result is for p = ^ : 



lim / uP{x,u,t)du 

.^oo J 



L 

n 



where K is the modified Bessel function of the second kind P^. Thus, 
uit,x) = ^-VT^ex\x\K{l + ^^7^) 

Calculations using asymptotic expansion of the K functions for e ^ show that 

u{t,x) = -X + 0(1), e^O. 

e 

Analogously result we get in the case of any rational p < I. Namely, if we assume 
P~m^' mi,TO2 G we obtain 



u{t,x) ^ ^x + o(^^^ , e^O. (4.7) 
In particular, ioi p — to G N 

u{t,x) = —X +0 \£ e — > 0. 

Then, we can consider p as a parameter and notice that at any fixed a; 7^ 0, e € 
(0, 1] the functions P{t{e),x, u;p) and uP{t{e), x, u;p) are continuous on the set u G 
f^L, p e [0, 1]. Moreover, ^ lim / P(t(e), X , u; p) du and J- lim J u P{t{e), x,u;p) du 

are uniformly bounded for p € [0, 1 — 5], (5 is a positive arbitrary small constant 
(see (4.6)). Thus, by a standard reasoning we conclude that the ratio u{t,x) is 
continuous in p, p S [0, 1 — (5]. Since there exist a rational sequence convergent to 
every real p, the continuity implies that (4.7) holds for all p E [0,1 — S]. Over the 
arbitrariness of S we can conclude that (4.7) is true for all p e [0, 1). 

Now we fix e > and expand the integrals in (2.2) near x — 0. We have for 
/? = 0, any fc = 1, n, using (2.9): 

Uk P {t{e) , X , u) du = 
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a2(e-l)(27ra2(e-l)a)"/^ " 

\aL\ 




= Cfe(e) Xk, 

for some functions Ck{s). Moreover, 



/ P(,W,x..)<i. = 



^ exp ] \ur('-^^-U\u\+o{\x\) 



2b 

All integrals converge for p > 0, p ^ 1 as L — > oo. Further, 

n 

due to the equality of space directions Cfe(e) = ^ ^(e), therefore 

= ja; + o(|a;|), |a;| ^ 0. (4.9) 

In the limit case p = I, where both integrals in the ratio (2.2) diverge, we can apply 
de L'Hopital's rule as above (in (4.9) to obtain (4.6)). 
Thus, the theorem is proved. 

□ 

Remark 4.1. For (3 > 0, p> 1 we get the expansion (4-3) as e — > again, hut in 
this case 

C = -i.n^ f2 r((„.2,(p 

+ aj J2ppj T(n{p-l)/2p) 



/3 a \ + 



forp e [0, 1) we get qualitatively the same result as in (4-V- 
Analogously we have for any p > 0, e G {0, 1] 

u{t,x) = — ^-g X + o(|a;|), |a;| 0. 



Remark 4.2. Let us compute the conditional variance v{t,x). For p > 1 + ^ we 
can expand the probability density P{t, u, x) in e near e = and obtain that for any 
x^O 

v{t{s),x) = F{n,p, a, a) \x\ i (l - + + o{e), 
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with 

/ „(p-l)-2 \ 

Thus, at any fixed x ^0 the variance tends to some finite value as £ — > 0. 

Further, for any fixed e > we get the following asymptotic expansion near 
x = 0: 

p/ n(p-l)-2 \ 2 \^ 

rdn) (- 4a^-e) ) +''<W)- 

therefore the variance tends to zero as x ^ and e — > 0. 

We cannot write an explicit formula for the variance for allp < 1 + ^- However, 
for p < 1 in the expression for the second moment both integrals in the ratio diverge 
as L ^ oo, and we can use de L'Hopital's rule once more to show that the variance 
tends to infinity as £ — > 0. In particular, for p = m e N, 



v{t,x)=o(e-^^ , 



0. 



5. Asymptotic behavior for a Gaussian initial distribution 



Let us set / = exp(— x ), fc > 0, and study the asymptotics of (2.2), 

(2.9) near = 0. We consider only the case (3 — and n = 1. As we have seen in 
Sec. 3, this type of initial density distribution eliminates the unbounded gradient 
growth as £ ^ in the case of a constant a. 

To get a qualitative result we firstly find an asymptotic expansion near a; = 
for both integrals in(2.2), and for their ratio in e at the point £ = {t = T). The 
computations made in the environment Maple show that for p < 1 

u = -Csx + o{e)Ci{x) + o{\x\)C2ie), 

and for p > 1 

U = -^X + Cx + 0{£)Ci{x) + 0{\x\)C2{£), 

with a positive constant C depending only on a, a, k and bounded functions Ci 
and C2. 

Numerical computations show that if p < 1 the function u{t, x) is very close to 
the respective function for p = 0. In particular, the singularity at the origin is also 
eliminated. Mp> 1, the picture is similar to the case of uniform initial distribution 

m- 

6. "Observable" density and "induced" velocity 
Introduce the function 



p{t, x) = lim / P{t, X, u) du. 

L^oo J 

If cr = 0, in (2.1), then p{t,x) solves the continuity equation 

1^ + div(wp) = 0, (6.1) 

and therefore corresponds to the distribution of density for the particles provided 
initial data p{0,x) = f{x) > are given. Here we do not require that f{x) is a 
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probability density. Thus, p{t, x) coincides with the real density of particles that 
we denote by p{t, x). 

For a the function pit, x) is not a solution to (6.1). Nevertheless an observer 
who does not know the random distribution of the particles positions considers 
p{t, x) as the density for the particles. We will call p{t, x) in this case an observable 
density. It differs from p{t, x). For example, it is easy to compute that for f{x) = C 
and (3 — G from the system (1.1), (6.1) one gets p{t,x) = = Ce~". The 

observable density p{t,x), as follows from (4.8), behaves in a quite different way. 
The only common feature is the asymptotics 0{e~") as \x\ 0. Nevertheless, given 
a density p, one can find from (6.1) the velocity u{t,x). Let us consider the case 
n= 1. If we assume 

lim p{t, x)u{t, x) =0, t€R+ 
(this implies the momentum conservation), we get 

X 

I Pt{t,x)dx 

u{t,x) = . 

p{t,x) 

Remark that p{t, x) is positive for f{x) > 0. Let us introduce the vector-function 
v{t, x) according to the formula 

X 

J p^{t,x)dx 

■ ^'-'^ 

It is natural to call v{t,x) the "induced" velocity. As follows from (6.2), (2.1), we 
have 

v{t,x) = u{t,x) + vi{t,x), (6.3) 

where 

2 / \u\^PP^it,x,u)du 

vi{t,x) := 



2 / P{t,x,u)du 

R 

One can see that if P{t, x, u) is given by (2.9) in the case f{x) = const the second 
term vi in (6.3) vanishes at a; = 0. Moreover, computations show that as |a;| — > 

vi{t,x) = C{t)sx\xf + o(|a;|^), 

with a function C{t) which stays bounded as e ^ 0. Thus, in the vicinity of the 
point of the singularity formation the induced velocity is close to u. 

7. Conclusion 

We have shown that in the model (1.8), generally speaking, the expectation of the 
velocity given position differs drastically from the velocity in the Burgers equation 
in the deterministic coordinate space. Under the observation in the random position 
of particles the blow up phenomena can be lost. In this context the uniform initial 
distribution of the particles position seems most interesting. Here the threshold 
effect arises. Namely, if the exponent p is less or equal then 1, then the expectation 
u follows the real behavior of the velocity rather well. However, ii p > 1, the 
function u tends to zero whereas the real velocity tends to infinity at any point 
outside the origin a; = as the time tends to the critical value. 
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The question on a PDE obeyed by the function u arises naturally, due to existing 
formalisms to represent solutions of PDE as the expected value of functionals of 
stochastic processes (see e.g. P1],|TS], [TB] and references therein). In TT one can 
find a recent result in the field close to our study, the stochastic formulation of the 
viscous Burgers equation. Namely, it was shown that if the pair {u,X) (where X 
is the flow map such that uo X is constant in time) solve the stochastic system 

dX = u dt + V2h' dW, u = 'E[uo o X^^], X{xo,0) — xq, v = const > 0, 

then under certain additional conditions u satisfies the viscous Burgers equation 

ut + (u, V)u = i^Au. (7.1) 

The function u also solve (7.1) for p — 0, however only for the special initial data 
and uniform initial distribution of particles. For example, one can check that u 
defined as (2.2), with P{t,x,u) governed by (2.1) does not satisfies the viscous 
Burgers equation even for p — and linear initial data for the Gaussian initial 
distribution (see (3.5)). Nevertheless, the comparison of u and the solution to the 
Burgers equation with a specific viscosity term is a very interesting open question. 

At last we would like to discuss possible applications of the results of this paper. 
First of all the the stochastic systems like (1.8) with velocity field of the fluid, with 
possibly random component with prescribed statistics can be applied to model 
turbulent or other disordered fluctuations. For example, the problem can be in 
describing some desired statistics of the trajectory X{t) of a tracer particle released 
initially from some point xq and subsequently transported jointly by the flow u{x, t) 
and molecular diffusivity |18j. 

Other applications can be found in the theoretical financial mathematics to 
model the earnings yield of a risky asset {X{t)) with the trend and volatility de- 
pending on a certain macroeconomic factor U{t) (e.g. the spot interest rate) in 
the spirit of [19]. Let us notice in this context that for p = in several cases it 
is possible to solve the Fokker-Planck equation for more general class of stochastic 
DE, describing the economic quantities in a more agreed way, namely 

m+n 

dXk{t) = (afe + AmU^H)) dt +Y1 diW^)t, k = 1, .., n; 

1=1 

m+n 

dU,{t) = {h + B,k Uk{t)) dt+Y, d{Wk)t, i - 1, .., m; 

fc=i 

Xk{0)^xk, U^{0)=u,, t>0, 

where W{t) is a M'"+" valued standard Brownian motion process, the market pa- 
rameters Ofc, Aifc, (Jifc, Bifc, are constant matrices of appropriate dimensions. 
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